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ABSTRACT 



A TOPAZ-II Thermionic Fuel Element (TFE) test stand is under study at the 
TOPAZ International Program. Extensive cross-sectional drawings of this system are 
presented along with detailed descriptions of its construction A thermal analysis of the 
test stand is developed using a thermal resistance network method and data collected on 
the test stand. Critical resistances in the network include the 0.5 mm regulated helium 
gap, cooling water channel, and cesium interelectrode gap. The thermal analysis indicates 
that for an input power to the active zone of 3003 W, and a regulated helium gap pressure 
of 10 torr, the thermal conductivity of the regulated gap helium ranges from 0.073-0.087 
W/mK. This is much less than the thermal conductivity calculated as a function of 
temperature only of 0.276 W/mK. For a water mass flow rate of 0.04 kg/s, the convective 
heat transfer coefficient ranges from 528-1012 W/m : K based on correlations for fully 
developed laminar and turbulent flow. Our model may be used in conjunction with a 
thermionic emission model in order to more completely evaluate the operating 
characteristics of the test stand, and to better correlate results with actual operating 
conditions in the TOPAZ-II reactor core. 
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I. INTRODUCTION 



The TOPAZ-II Thermionic Fuel Element (TFE) test stand was originally built and 
tested at the Scientific Industrial Association "LUTCH" in Podolsk, Russia In December, 
1992, the stand was reassembled at the New Mexico Engineering Research Institute 
(NMERI) in Albuquerque, New Mexico. The TFE test stand was designed to conduct 
single TFE non-nuclear testing by replacing the nuclear fuel with a tungsten heating 
element (TISA heater). The test stand provides a means to test a single TFE so that its 
performance can be better observed without the interference and complexity of the 
integrated TOPAZ-II system. This removes the direct coupling of the TOPAZ-II in-core 
TFE with the coolant loop, radiator and the other TFEs electrically connected in series. 
Acceptance and demonstration testing were conducted throughout 1993 (Luchau et al., 
1993). 



In order to better understand and predict the behavior of a TFE in the test stand, a 
detailed thermal analysis based on a network method needs to be developed to account for 
the thermal losses in the stand. Because of the large amounts of Russian technical data 
and engineering drawings associated with the TOPAZ-II TFE test stand, and due to the 
extreme difficulty in translating Russian technical information, the implementation of this 
thermal model required the creation of detailed and extensive cross-sectional drawings for 
this thesis. Mathcad Plus 5.0 (MathSoff, Inc., 1994) was the primary technical calculation 
software used to solve the thermal analysis model. 

Since the test stand provides a means for testing a single TFE without the support 
of the supporting TOPAZ-II reactor systems, there are several differences in the 
construction and operation of the test stand from the TOPAZ-II reactor core. In the test 
stand, water cooling is used to remove the excess heat instead of the liquid metal cooling 



1 



that is used in the reactor core. The cooling water flows through the test stand water 
jacket along a torturous path to aid in its heat removal ability. The heat transfer 
coefficient for the cooling water is directly dependent on the flow cross-sectional area and 
the total heat transfer area, both of which are challenging to determine. Thermocouples 
internal to the test stand provide a axial temperature distribution along the TFE working 
section. No internal temperature measurements are made in the TOPAZ-II reactor core 
TFEs. Heaters are also located internally in the test stand to provide a means of degassing 
the test stand components prior to operation. A 0.5 mm gap in the test stand just inside of 
the cooling water jacket is filled with low pressure helium to provide a means of 
controlling the heat transfer through the test stand. This gap is used to help the test stand 
better simulate the operating conditions experienced by a TFE in the reactor core. Due to 
the low pressure band that the gap is adjusted over (from 1 to 10 torr), the pressure 
dependence of the rarefied helium's thermal conductivity is difficult to define. 

Experimental data for this work were collected in March 1994, with follow-on 
experimental data taken in August 1994 at NMERI. Along with aiding in the development 
of the thermal model for the test stand, these data are used to investigate several of the 
critical resistances from the thermal network. A good understanding of the test stand 
thermal characteristics is essential in evaluating the performance of a single TFE operating 
in the test stand, and correlating that performance to a TFE operating in the TOPAZ-II 
reactor core. The experimentally measured input power to the TFE working section along 
with the collector sleeve temperature information is also used as an input to the TOPAZ-II 
Thermionic Transient Analysis Model (TITAM). TITAM is a computer simulation 
program developed by the University of New Mexico to look at both steady-state and 
transient operating characteristics of the TOPAZ-II reactor system (Paramonov -and 
El-Genk, 1994). 
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The goals of this work include the creation of detailed engineering drawings for 
the test stand, investigation of test stand material properties, development of a thermal 
model for the test stand, and benchmarking the model with the experimental data 
collected. Additionally, experimental data are used to empirically fix a range for several of 
the thermal properties that cannot be easily determined 

LT J. R. Venable is writing a companion thesis that will concentrate on the 
thermionic theory behind the data collected for this thesis. Venable's thesis will be 
available after March of 1995 and a joint paper will be presented ("Operational Testing 
and Thermal Modeling of a TOPAZ-II Single-Cell Thermionic Fuel Element Test Stand") 
at the Space Nuclear Power and Propulsion Symposium held in Albuquerque, New 
Mexico, on January 8-12, 1995. 
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II. BACKGROUND 



The continued exploration and exploitation of space will require increasing 
amounts of power as technology continues to advance. Supporting high definition 
television broadcasts and the rapidly growing world-wide communication needs, moving 
large level research and industrial opportunities to Earth orbit or into the cislunar region, 
revisiting and colonizing the Moon, and pushing the envelope of civilization to Mars and 
beyond are future goals which will demand large electrical power sources that are both 
compact and have long lifetimes. These larger power sources will need to survive harsh 
environments, including radiation belts, the vacuum of deep space, and hostile 
atmospheres of other planets and natural satellites in the solar system. Nuclear power can 
provide the compact power sources necessary to meet these future needs. 

A. UNITED STATES SPACE NUCLEAR POWER PROGRAMS 

Since the mid-1950s, the United States has been interested in the use of nuclear 
power in space applications. From radioisotope power systems, nuclear rockets, and the 
use of nuclear reactors for power production, interest in the space nuclear power area was 
motivated by a favorable political environment and public support until the early 1970s. 

The Transit 4 A navigational satellite was launched in June of 1961 with a 
SNAP-3B Radioisotope Thermoelectric Generator (RTG) to provide power. Since that 
first successful launch of a nuclear power source, a total of 24 launches using RTGs have 
followed including the Pioneer, Viking, and Voyager missions. Although the RTG has 
proved to be very useful for many applications, the ability to support a very ambitious 
space program required the development of a nuclear reactor based on the principle of 
atomic fission for both propulsion purposes and to provided electrical power. 
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Project Rover was initiated at Los Alamos National Laboratory in 1955 and was 
the beginning of approximately 18 years of research and development into nuclear rockets. 
Although the nuclear rocket program was terminated in late 1972, it was the start of the 
search for a high temperature nuclear reactor The research and demonstration test of the 
Nuclear Engine for Rocket Vehicle Application (NERVA) program proved that a nuclear 
reactor could provide a high power density rocket with high thrust and specific impulse. 
(Angelo and Buden, 1985) 

Active development of a nuclear reactor to provide electric power began in 1957 
with the Systems for Nuclear Auxiliary Power (SNAP) program. These nuclear reactors 
required a power generation method to convert the heat created by atomic fission into 
electrical power. The initial concept was SNAP-2, a thermal uranium zirconium-hydride 
reactor that utilized a Rankine cycle with mercury as the working fluid to convert the 
nuclear heat to electricity. This was the precursor to two separate concepts, the SNAP-8 
and the SNAP-10A, both of which were initiated in 1960. The SNAP-8 reactor utilized a 
mercury Rankine cycle along with a larger version of the SNAP-2 reactor system to 
support anticipated higher power demands. (Angelo and Buden, 1985) 

The SNAP-10A was also based on the SNAP-2 reactor system, but instead of 
using the dynamic Rankine power conversion cycle, it used a thermoelectric power 
conversion system. The use of static power conversion provides a means of converting 
heat directly into electricity without the need of a secondary working fluid (such as 
mercury for a Rankine cycle or helium-xenon for a closed Brayton cycle) and the 
associated turbines and pumps/compressors. In April of 1965, the SNAPSHOT mission 
was successfully launched carrying a SNAP-1 0A reactor to approximately a 13 00 'km 
orbit. Although the SNAP-10A reactor shutdown only 43 days after becoming critical, it 
proved that a nuclear reactor could be successfully launched and operated on orbit. This 
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was the only attempt by the United States to launch a nuclear reactor into space. (Angelo 
and Buden, 1985) 

Another static power conversion method that was developed in parallel with the 
thermoelectric method was thermionic power conversion. Whereas thermoelectrics work 
on the basis of a thermally induced emf at the junction of two dissimilar semiconductors, 
thermionics works basically by "boiling off' electrons from a hot emitter surface and 
collecting the electrons on a relatively cooler collector surface across a small gap. The 
collector and emitter are then electrically connected together through an external load 
which completes the electrical circuit. The process of thermionics is discussed in more 
detail in the next chapter. Although thermionic devices were tested in the mid- to 
late- 1960s, no complete reactor systems using thermionics was actually built in the United 
States. (Angelo and Buden, 1985) 

Due to a changing political environment and shifting emphasis in the space 
program, all nuclear reactor development programs were canceled at the end of 1972. 
Although some research did continue after 1972, it wasn't until the early 1980s that 
increased interest was returned to space nuclear reactors. This renewed interest was 
mainly due to the need for a large power source to support the Strategic Defense Initiative 
(SDI). The advance space nuclear power program was reinitiated as the Space Power 
Advanced Reactor (SPAR) system later renamed SP-100. Several reactor designs and 
power conversion systems were initially considered for SP-100. The primary concept that 
SP-100 developed utilized a thermoelectric power conversion system. (Angelo and 
Buden, 1985) Research on SP-100 has continued until 1994 when once again the political 
environment has changed and funding for the program is being terminated. 
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B. 



FORMER SOVIET UNION SPACE NUCLEAR POWER PROGRAMS 



The former Soviet Union, like the United States, was very active in the research 
and development of space nuclear power. From the early 1960s to 1989, a total of 38 
nuclear reactor systems were launched by the former Soviet Union. The majority of the 
nuclear reactors launched were for the Radar Ocean Reconnaissance SATellites 
(RORSATs) and utilized thermoelectric energy conversion. The other two nuclear 
reactors launched were COSMOS-1818 and COSMOS- 1867. Both COSMOS satellites 
were test platforms for the TOPAZ-I reactor system utilizing multi-cell thermionic energy 
conversion. (Polansky et al., 1994) 

In the early 1960s, the former Soviet Union began development of single-cell and 
multi-cell TFEs. In 1967 the two concepts were sent for system development to two 
separate technical institutes. The single-cell TFE was developed by the Central Design 
Bureau of Machine Building (CDBMB) in St. Petersburg under the project name of 
ENISY. The multi-cell TFE was developed by the Scientific and Production Cooperation 
"Krasnaya Zvezda" in Moscow under the project name of TOPAZ, the Russian acronym 
for Thermionic Experiment with Conversion in Active Zone. Both concepts were 
developed cooperatively between the two institutes up until the early 1970s. This lead to 
many similarities between the two reactor designs. Cooperative efforts ended in the early 
1970s due to differences between the two institutes and the programs developed 
independently. (Voss, 1994) 

In 1987 and 1988, the multi-cell TFE reactor TOPAZ was successfully launched 
twice on COSMOS-1818 and COSMOS- 1867. The first mission operated for 142 days 
and the second mission operated for 342 days before the planned depletion of the cesium 
supplies on board terminated the tests. Both test were extremely successful and confirmed 
that the on-orbit data correlated well to ground test results. (Paramonov and El-Genk, 
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1994) Although there were 26 single-cell TFE ENISY reactor systems manufactured and 
19 systems tested between 1970 and 1989, none of the single-cell TFE technology were 
ever launched and tested on orbit. (Schmidt et al., 1994) 

C. TOPAZ INTERNATIONAL PROGRAM 

At the end of the 1980s, the Russians approached a private U S. company 
(International Scientific Products, San Jose, CA) about transferring Russian thermionic 
technology through the purchase of Russian space reactor hardware. During the 
negotiations, the U S. company inadvertently renamed ENISY to TOPAZ-II and TOPAZ 
to TOPAZ-I (Voss, 1994). In May of 1992, the U.S. Government purchased two 
complete TOPAZ-II reactor systems and associated test and support hardware through 
the purchasing agent International Scientific Products. Funding was provided by the 
Strategic Defense Initiative Organization (SDIO) to investigate the possible use of the 
TOPAZ-II reactor system in support of the Nuclear Electric Propulsion (NEP) Space Test 
Mission. This project has since been named the Topaz International Program (TIP). 
(Wyant, 1994) 

The TEP facility is located at NMERI which is owned and operated by the 
University of New Mexico in Albuquerque, NM. Funding is still being provided by the 
former SDIO (now the Ballistic Missile Defense Organization, BMDO) and administered 
by the U.S. Air Force's Phillips Laboratory. The New Mexico alliance for Thermionics 
which includes the University of New Mexico, Sandia National Laboratories and Los 
Alamos National Laboratory provides research support for the TIP program. The 
international effort includes scientist and engineers from the Russia, England, and Fraince. 
Although the flight program for TIP has been canceled, in March of 1994 four additional 
TOPAZ-II units were purchased for continuing research into thermionics and reactor 
operation. (Wyant, 1994) 
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III. THERMIONIC ENERGY CONVERSION 



A. BASIC THEORY 



The thermionic process can be considered to be a "heat engine" which utilizes 
electrons as its working fluid to convert heat into electricity. The thermionic converter 
consists of a emitter, collector, and a small gap between the two. The collector is 
electrically connected to the emitter through an external load to complete the electrical 
circuit. Figure 3-1 shows a simplified flat-plate thermionic converter (Rasor, 1991). 
Electrons are "boiled" off the hot emitter and travel across the gap and are "condensed" on 




p = q — Q 

LOAD IN LOSS 

Figure 3-1 . Simplified Flat-plate Thermionic. 
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the relatively cooler collector The potential difference between the emitter and collector 
creates a driving force which causes a current to flow through the external load. (Angelo 
and Buden, 1985) 

Since the thermionic device acts as a "heat engine," the maximum efficiency of the 
thermionic cycle can be characterized as the cycles Carnot efficiency. Some of the 
inefficiencies of the thermionic cycle include: (1) radiation heat transfer from the emitter 
to the collector; (2) negative space charge in the gap that interferes with the electron flow 
from the emitter to the collector; and (3) thermal losses. These inefficiencies prevent the 
thermionic cycle from approaching Carnot cycle efficiency. (Angelo and Buden, 1985) 

Radiation heat transfer from the emitter to the collector cannot be mitigated. 
Higher collector temperatures would not only reduced the amount of radiation heat 
transfer but would also allow for a much smaller radiator since the radiators surface area is 
inversely proportional to the radiators operating temperature to the fourth power. 
However, as the collector temperature increases, the number of electrons being "boiled" 
off the collectors surface, called electron back emission, also increases. This in effect 
causes the overall cycle efficiency to decrease dramatically. A design trade-off must be 
made to understand the optimum collector operating temperature. (Angelo and Buden, 
1985) 



The negative space charge in the gap can be decreased by either making the gap 
size smaller or by injecting a low pressure gas into the gap. A vacuum diode is a 
thermionic device that relies on a very small gap between the electrodes to reduce the 
current inhibiting effects of the electric field produced by the negative space charge. 

Theoretical spacing needed to produce useful current density approaches the associated 
Debye length (~ 1 pm) which is for a practical purpose very difficult to manufacture and 
maintain. This very narrow gap can be widened somewhat (up to at most 0.5 mm) by 

injecting a gas plasma into the gap. This type of thermionic device is known as a plasma 
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diode. The primary purpose of this plasma is to facilitate to transfer of electrons through 
the gap while minimizing the potential drop and current attenuation due to the presence of 
the plasma. (Rasor, 1991) 

The most common gas that is used in this type of diode is cesium which has the 
lowest ionization potential of any chemical element (Angelo and Buden, 1985). A 
measure of the amount of energy necessary for a electron to escape through the potential 
barrier associated with the surface of a metal is called the work function and is given in 
units of electron-volts (eV) (Angrist, 1976). Not only does the cesium plasma in the gap 
decrease the space charge, it also may have the effect of forming a monolayer of cesium 
on both the emitter and collector causing a decrease in their associated work functions. 
This decrease in the metals work function is directly related to the pressure of the cesium 
in the gap and the electrode temperature. Cesium pressure can be controlled by adjusting 
the temperature of the liquid cesium reservoir which in effect controls the vapor pressure 
of the evaporating cesium. (Angelo and Buden, 1985) 

B. TOPAZ-II NUCLEAR REACTOR 

The TOPAZ-II nuclear reactor consists of 37 individual TFEs incorporated in a 
cylindrical core arrangement. Figure 3-2 shows a cross-sectional view of the TOPAZ-II 
reactor. Each of 37 TFEs consist of a UO ; fuel assemble, a single-cell cylindrical 
thermionic device, and a sodium-potassium (NaK) liquid metal cooling channel. Electrical 
power is supplied to the loads by 34 of the TFEs connected in series. These TFEs are 
designed to produce 5-6 kW e at 27 VDC for up to 3 years. The other three TFEs are 
connected in parallel to provide power to the electromagnetic pump. This pump is used to 
generate the pumping power necessary to circulate the NaK. The overall length of the 
reactor core is 37.5 cm with an overall width of 26.0 cm. (Paramonov and El-Genk, 
1994) 
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Figure 3-2. Cross-sectional View of the TOPAZ-II. 

The TFEs are inserted into axial channels in a zirconium-hydride (ZrH) moderator. 
Twelve reflector drums surround the core, three are safety drums which are turned out 
prior to initial criticality when on orbit. The other nine reflector drums are used to control 
the nuclear reaction during operation and are called control drums. All of the reflector 
drums are made out of beryllium clad in stainless steel with a small section of neutron 
absorbing material made of boron-carbide (B 4 C). A beryllium reflector surrounds the 
remainder of the core as shown in Figure 3-2. (Schmidt et al., 1994) 

Excess heat is removed from the TFEs using a eutectic NaK coolant pumped 
through stainless steel piping by the electromagnetic pump. The hot coolant flows 
through a conical radiator where the waste heat is rejected to space. During normal 
operation the thermionic emitter temperature ranges between 1800 K to 2100 K and the 
NaK coolant maintains the collector temperature between 743 K to 843 K. (Schmidt et 
al., 1994) 
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c. 



TOPAZ-II THERMIONIC FUEL ELEMENT 



The TOPAZ-II TFE consists of a single-cell cylindrical thermionic converter. A 
single-cell converter utilizes the entire length of the reactor core as a single device. This is 
in contrast to the TOPAZ-I TFE which is made up of multi-cell thermionic converters. A 
simple comparison of the electrical circuit for both a single-cell and multi-cell TFE are 
shown in Figure 3-3. The advantages of using multi-cell converters are that they produce a 
higher power density and higher conversion efficiency. This is mainly due to less ohmic 
losses in the multi-cell because of the shorter thermionic elements. Another advantage of 
the multi-cell design is that each of the cells are connected in series resulting in a TFE 
output voltage in the volt range. This reduces the amount of electrical losses in the leads 
and makes the power distribution system more efficient. (Paramonov and El-Genk, 1994) 



/Electron \ /Electron \ 

'Current ' 'Current ' 





MULTI-CELL TFE SINGLE — CELL TFE 

Figure 3-3. Simplified Electrical Schematics for Multi-Cell and Single-Cell TFEs. 
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The two major disadvantages of using a multi-cell TFE are that manufacturing is 
significantly more complicated and the design does not provide a means of removing the 
nuclear fuel and performing non-nuclear testing. Although the single-cell TFE design has 
a much lower efficiency, it does provide a less complicated converter and a means of 
performing non-nuclear testing. This testing is accomplished by replacing the UO ; fuel 
assembly with a tungsten heater (TISA heater). (Paramonov and El-Genk, 1994) The 
TISA heaters provide a temperature distribution similar to that created using nuclear 
heating. Detailed information on the use of electrical heaters to simulate nuclear heating 
can be found in El-Genk and Xue, 1 994. 

The interelectrode gap between the emitter and collector is filled with cesium 
vapor during normal operation. The cesium is supplied by a single liquid cesium reservoir 
that is heated by one of the NaK coolant return lines. The cesium in the gap is constantly 
vented to space at a rate of approximately 0.5 g/day. This is one of the limiting factors in 
the reactor lifetime. (Schmidt et al., 1994) During normal operation, the cesium pressure 
in the gap is maintained around 1. 5-2.0 torr (Paramonov and El-Genk, 1994). 

Axial and radial temperature distributions in the reactor core cause performance 
differences in the 37 TFEs. The exact performance characteristic of each individual TFE 
is difficult to determine due to the following: (1) the thermal coupling of each TFE with 
the TFEs surrounding it, or lack of TFEs; (2) the electrical coupling of each TFE with the 
other TFEs, either in series for the 34 load power producing TFEs or in parallel for the 3 
pump power producing TFEs; (3) all TFEs share a common heat removal system; and (4) 
all TFEs share a common cesium system. To understand the individual thermionic 
performance of a TFE, the TOPAZ-II TFE test stand must be used. This test stand allows 
detailed analysis of a single TFE under unique operating characteristics. 
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IV. TEST STAND DESCRIPTION 



The TOPAZ-II TFE test stand is used to perform non-nuclear testing of a single 
TOPAZ-II TFE. This allows the TFE performance to be determined without the 
interference of the rest of the TOPAZ-II reactor systems and without any interfering 
effects caused by the TFEs relative location in the core (see Figure 3-2). Interference by 
the reactor systems include the direct thermal coupling of the TFE to the NaK coolant, 
coolant pump, and radiator, and the electrical coupling of the TFE to the other TFEs 
connected in series in the reactor core. The test stand also allows varying certain 
parameters which could not be varied on the entire reactor system to investigate their 
effect on TFE performance. These parameters include varying cesium vapor pressure in 
the interelectrode gap, the removal of all cesium vapor from the interelectrode gap, 
changing the collector sleeve temperature (directly effecting the collectors temperature), 
and operating the TFE at its operational limits to determine break-down characteristics. 

The test stand contains the test apparatus within a vacuum chamber along with the 
associated vacuum, cesium, helium, and cooling water systems external to the vacuum 
chamber. As with almost all Russian structural components, the vacuum chamber is 
constructed mainly of Russian stainless steel type 12X18H10T. Phillips Laboratory has 
determined that this type of Russian stainless steel is closely approximated by the 
American Iron and Steel Institute (AISI) standard 321. The vacuum system is used to 
initially evacuate the vacuum chamber, the helium cavities, the cesium vapor system, and 
the TISA heater channel. The test stand start-up procedure requires an involved 
evacuation and out-gassing procedure. Out-gassing is performed by numerous heaters 
throughout the test stand. 
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The cesium system is used to supply cesium vapor to the interelectrode gap and to 
maintain the cesium vapor pressure. The system also provides a means to purify the 
cesium inventory. When the test stand is operating, cesium vapor is provided to the 
interelectrode gap by a cesium reservoir that is inside the vacuum chamber. This isolates 
the internal cesium system from the rest of the system outside of the vacuum chamber. 
Cesium vapor pressure in the interelectrode gap is regulated and varied by adjusting the 
inner cesium reservoir saturation temperature. An external cesium reservoir is provided to 
evacuate and fill the cesium system when the test stand is not in operation. (Luchau et al., 
1993) 



The helium system provides helium to the cavities and gaps inside the test stand, a 
means to purify the helium supply and a means to evacuated the helium from the system. 
Helium is used to fill in the gaps in the test apparatus to facilitate and control the heat 
transfer process. Collector sleeve temperatures can be controlled by varying the pressure 
in the regulate helium gap. The mechanism by which the thermal conductivity of helium is 
affected by pressure is discussed in detailed in the next chapter. 

The water cooling system picks up the excess heat throughout the test stand. The 
components cooled by the water system include electric drives for various valves, vacuum 
pumps, cesium reservoirs, transistor load bank and the test apparatus. The volumetric 
flowrate of the cooling water flow to the test apparatus is measured using a turbine type 
flowmeter in the inlet line. The inlet and outlet cooling water temperature is measured by 
type-K, surface contact thermocouples on the inlet and outlet stainless steel piping. Two 
thermocouples are used on the outlet piping, one on the outlet pipe that is in thermal 
contact with the vacuum chamber and the other on the outlet pipe that is not in thefmal 
contact with the vacuum chamber. The experimental data show that both of these 
temperature measurements are well within the thermocouple tolerance of each other. The 
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flowrate and temperature measurement accuracies are critical to the thermal modeling of 
the test stand. 

The stand utilizes an automatic control system, an automated data acquisition 
system, and a variable electrical load bank. The entire control and data acquisition system 
was completely upgraded from the original Russian system to a United States technology 
based system. This upgrade occurred between the first set of tests in March of 1994 and 
the second set in August of 1994; however, the testing parameters were maintained the 
same during both data sets. The load bank is composed of both a transistive and a 
resistive load bank. The load bank can be controlled by a micro-computer to provide 
current-voltage (I-V) sweeps for the TFE. These sweeps are important in determining the 
actual operating characteristics of the TFE. (Luchau et al., 1993) 

A diagram of the test apparatus is shown in Figure 4-1 (based on Luchau et al., 
1993). The upper and lower helium chambers provide helium to the unregulated gaps in 
the test section and to the spaces in the TISA heater. The cesium for the inter-electrode 
gap enters the test section through the lower helium chamber which is heated during 
testing to prevent "cold-spots" in the cesium line. Since the lower helium chamber is 
heated and the upper helium chamber is not heated, heat losses out of the top and bottom 
of the test apparatus are different. This difference in the end losses will affect the test 
stand's thermal model. The upper and lower helium chambers also contain the electrical 
connections to both the emitter and collector. Electrical connections to the TISA heater 
are made in a chamber above the upper helium chamber. 

A lengthwise cross-sectional view of the test stand's working section is shovfti in 
Figure 4-2. This figure shows the relative proportions of the various test section 
concentric cylinders. The overall length of the thermionic working section is 375.0 mm. 
The length of the TISA heater is 300.0 mm and the length of the cooling water jacket is 
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Figure 4-1 . Cross-sectional View of the TOPAZ-II TFE Test Stand. 
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Figure 4-2. Cross-sectional View of a TOPAZ-II TFE in the Test Stand. 

(all deminsions in mm) 



approximately 326.0 mm. The extremely small dimensions and close tolerances gives an 
indication of the difficulty associated with the manufacturing of the test section. A radial 
cut-away view of the complete test stand is shown in Figure 4-3. An enlarged view of the 
portion of the test stand from the collector sleeve inward is shown in Figure 4-4. The 
overall diameter of the entire test section is only 87.0 mm giving an L/D of 4.31 for the 
working section. 

The TFE and its associated TISA heater are fit into the collector sleeve tube of the 
test stand as shown in Figure 4-4. The TISA heater is constructed of a special tungsten 
heating element and flexible lobes that act as thermocompensators. The maximum supply 
voltage is 29 VAC with a maximum current of 170 A. The TISA heater can supply a 
maximum heating power to the working section of 4500 W (Wold et al., 1993). The inner 
conductor of the TISA heater has a diameter of 6.5 mm and the outer conductor has a 
thickness of 0.4 mm and an outer diameter of 7.0 mm. The length of the heater is 300.0 
mm which makes the heater 37.5 mm shorter at each end than the actual thermionic 
working section. This shorter length simulates the end effects associated with nuclear 
heating. Detailed discussions of the correspondance of TISA heating to nuclear fuel 
heating are presented in El-Genk and Xue (1993). 

The thermionic device consists of concentric emitter and collector cylinders 
separated by an interelectrode gap and is 375.0 mm in length. The emitter cylinder is 1.15 
mm thick with an outer diameter of 19.6 mm. The emitter is made of a monocrystalline 
molybdenum - 3% niobium alloy with its exterior surface coated with a 0.1 mm thick 
monocrystalline tungsten layer. The 0.5 mm interelectrode gap is filled with cesium vapor 
during normal operation. The interelectrode gap is maintained by the use of 
scandium-oxide (Sc 2 0 3 ) spacers which prevent shorting between the emitter and collector. 
A total of six spacer sets are located axially along the TFE working section. These 
spacers are not shown in any of the figures. The collector cylinder is 1.4 mm thick with an 
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Figure 4-3. Radial Cut-away View of a TOPAZ-II TFE in the Test Stand 



Figure 4-4. Radial Cut-away Detail of a TOPAZ-H TFE in the Test Stand. 
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Stainless Steel 




inner diameter of 20.6 mm. The collector is made of polycrystalline molybdenum with its 
outer surface coated by a 0.15 mm layer of alumina (A1 : 0 3 ). The alumina electrically 
insulates the TFE from the remainder of the test rig. The thermionic device and TISA 
heater used in the test stand are identical in construction to that used in the actual 
TOPAZ-II reactor. 

A 0.05 mm gap separates the alumina insulator from the collector sleeve tube. 
This gap allows the TFE to be fit into the test rig and is filled with helium during normal 
operation. The helium is provided by the upper and lower helium chambers (see Figure 
4-1). The helium pressure in this gap is unregulated and is normally in the 200 to 300 torr 
range. 



The collector sleeve tube is 3.0 mm thick with an outer diameter of 29.9 mm. It is 
constructed ofRussian 1X18H10T stainless steel. The collector sleeve contains 12 evenly 
spaced, 2.0 mm deep slots in its outer surface that run the length of the collector sleeve. 
Figure 4-5 shows the relative location of the thermocouple slots in the collector sleeve. 
These slots contain the only thermocouples located in the interior of the test rig. The 
thermocouples are positioned along the collector sleeve to provide an axial temperature 
distribution in the working section. The thermocouple wires are packed in alumina inside 
of stainless steel tubes which are laid inside the slots. These thermocouples are positioned 
in the test rig at the location where the NaK liquid metal cooling flows through the TFE in 
the actual TOPAZ-II reactor. This allows correlation of the test stand data to data 
obtained from the complete reactor testing. Another 0.05 mm helium gap separates the 
collector sleeve tube from the copper conductor. This small gap is for assembly purposes 
only. The helium in this gap is also supplied from the upper and lower helium chambers 
and is unregulated in pressure. 
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Figure 4-5. Thermocouple Slot Location in the Collector Sleeve. 
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The copper conductor is 15 0 mm thick with an outside diameter of 60.0 mm. It is 
made of 99.99% pure copper and provides a heat transfer medium outward to the cooling 
water. The outer surface contains five radially spaced slots that each contain four 
nichrome (NiCr) heating elements. These heaters are used for outgassing the rig prior to 
startup and are secured during normal operation The heater wires are insulated with 
alumina and sheathed in stainless steel. A 0.5 mm gap separates the copper conductor 
from the inner cooling water jacket. This gap is filled with helium supplied from a source 
external to the test rig. The pressure of the helium is regulated in the band of 1 to 10 torr 
during normal operation to control the collector sleeve temperature by varying the thermal 
resistance of the gap. The exact mechanism through which the helium thermal 
conductivity is affected by helium pressure is discussed later in the thesis. 

The cooling water jacket is made of Russian 1X18H10T stainless steel and has an 
inner wall thickness of 2.5 mm and an outer wall thickness of 1 mm. A coiled stainless 
steel tube between the jacket walls provides a spiral flow path for the cooling water. The 
water flows from the bottom to the top of the water jacket around the outside of the coil 
in the cooling "channel" depicted in Figure 4-2. The coils are spaced by 35.0 mm and 
spiral approximately 6.5 times from bottom to top. Cooling water temperature is 
measured on the inlet and outlet piping exterior to the vacuum chamber using surface 
thermocouples. Cooling water flow rate is measured in the inlet piping using a 
turbine-type flowmeter located exterior to the vacuum chamber. 
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V. TEST STAND THERMAL ANALYSIS 



As discussed in Chapter III, the thermionic cycle is a relatively inefficient 
conversion process. The TOPAZ-II TFE has a maximum efficiency of 5-6%; therefore, at 
least 94% of the heat generated by the TISA heater must be dissipated as a loss. Detailed 
thermal analysis of the test stand can provide a better understanding of how testing a 
single TFE in the test stand correlates to operating a TFE in the TOPAZ-II reactor core. 
The developed thermal model also yields valuable information about the limitations of the 
test stand for future testing and the validity of past and future test results. 

A. THERMAL MODEL INPUTS 

The amount of input power to the test stand is determined by measuring the input 
current to the TISA heater through an external shunt and by measuring the voltage drop 
across the input terminals of the TISA heater. Multiplying the current in amps by the 
voltage in volts gives the input power to the TISA heater in watts and is designated as 
P-hsa- This TISA input power must be corrected to calculate the amount of power actually 
delivered to the emitter of the TFE. Correction is required to take into account the lead 
and line losses associated with the heater. Although Russian documentation for 
calculating the correction factor was not available, the Russian technicians at TIP use a 
correction factor of 88%. Correcting the input power to the TISA heater by this factor 
provides the amount of actual power delivered to the "active zone" of the TFE. This 
active zone power is designated as P u . 

The output power produced by the TFE is determined by measuring the current 
through the external load bank and by measuring the voltage drop across the upper and 
lower emitter or collector terminal connections (see Figure 4-1). Since the voltage drop is 
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measured directly across the interior terminal connections of either the emitter or 
collector, no correction factor is required to determine the output voltage of the TFE. 
Multiplying the voltage in volts by the current in amps gives the TFE output power in 
watts and is designated as P out . 

The only temperatures that are actually measured internally in the test stand are the 
collector sleeve temperatures. These measurements occur 1.8 mm into the stainless steel 
collector sleeve from its inner cylindrical surface (see Figure 4-4). As discussed in 
Chapter IV, the thermocouples in the collector sleeve provide axial temperature 
measurements along the length of the working section. There are a total of 12 
thermocouples in the collector sleeve which are designated T56 through T67. During the 
data collection runs, thermocouple T64 was inoperative and was disregarded in the 
analysis. Thermocouples T59 and T60 and thermocouples T61 and T62 measure the 
temperature at the same location along the collector sleeve. Relative locations of the 
thermocouples with respect to the centerline of the working section are shown in Figure 
5-1. The lengths of the TISA heater, the thermionic working section, and the water jacket 
and the approximate locations of the cooling water inlet and outlet are also shown in 
Figure 5-1. 

In a one-dimensional model, it is necessary to take an average of the 1 1 operating 
collector sleeve thermocouples to obtain a value to use as an input temperature. Although 
many variables effect the magnitude of the collector sleeve axial temperature profile, the 
profile shape is generally the same under all conditions. Figure 5-2 shows a typical axial 
temperature profile at a TISA heater input power of 2400 W. The temperature data 
points are shown as diamonds with T56 at the -205 mm position and T67 at the 205' mm 
position. The dashed line represents a third-order polynomial curve fit of the data. 
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Figure 5-1 . Axial Location of the Collector Sleeve Thermocouples. 
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Figure 5-2. Typical Axial Temperature Distribution along the Collector Sleeve. 

Although the TISA heater provides a fairly uniform temperature profile into the emitter, 
several reasons inherent to the test stand design cause the distribution to take the shape as 
shown in Figure 5-2. An overall downward slope in the data is expected due to the 
cooling water flowing from the bottom to the top of the working section; however the end 
effects of the upper and lower helium chambers cause the data to not be linear at the two 
ends of the graph. A straight average of the 1 1 data points would bias the temperature 
input to the thermal model low due mainly to the end effects. This biasing can be 
minimized if the end data points (T56 and T67) are not used and an average is taken only 
of the remaining nine data points. For example, the average of the 1 1 data points in 
Figure 5-2 is 379.0 K. Excluding the end data points gives an average of 392. 1 K which is 
approximately 3.5% higher. For all calculations performed in this analysis, the average 
collector sleeve temperature was determined by excluding T56 and T67. A much more 
detailed two-dimensional model would be required to understand the error associated with 
this approximation. 
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The cooling water system provides the last of the measured test stand quantities 
necessary for the thermal analysis. Instrumentation of the cooling water system was 
required in support of this thesis since the information on the cooling water system was 
not previously deemed necessary. Surface thermocouples were placed on the cooling 
water inlet and outlet lines, exterior to the vacuum chamber. Approximately one meter of 
piping separates the actual thermocouple location from the inlet and outlet piping 
penetrations on the test stand. No corrected factor was used on these temperature 
measurements for any potential piping losses between the test section and the exterior of 
the vacuum chamber because the only method of heat transfer from the water piping 
would be by radiation heat transfer to the interior of the vacuum chamber. Since the 
maximum water temperature inside the piping is less than 40 °C over all the data runs, 
radiation heat transfer is not significant. The water piping is also in thermal contact with 
the vacuum chamber which provides another potential heat loss by conduction heat 
transfer to the vacuum chamber itself. Some uncertainty could have been removed if the 
thermocouples had been placed on the piping inside the vacuum chamber at the entrance 
and exit of the test section. This was not done for any of the data runs used in this 
analysis due to time constraints and the complexity of removing the vacuum chamber. An 
average bulk water temperature in the water jacket can be estimated by taking an average 
of the inlet and outlet water temperatures. 

A turbine-type flow meter was also installed on the inlet cooling water line prior to 
the data runs. The accuracy of the flow meter was physically verified to be within 2% of 
the indicated flow rate over the range used during the data runs. Initially the data runs 
were taken with the cooling water isolation values fully open which was the normal 
operating procedures for the test stand. However, with the valves fully open, the change 

in cooling water temperature between the inlet and the outlet was less than 5 °C. Since 
the accuracy of the surface thermocouples is only ± 1 °C, a larger temperature change is 
desired to reduce the percent error. Around halfway through the first set of data runs, the 

flow rate was throttled down to approximately 0.7 gpm. This gave a temperature change 

33 



between 10 to 15 °C over most of the data runs. All of the data collected during the 
second set of data runs were collected at a flow rate of approximately 0.7 gpm 

B. THERMAL MODEL 

Difficulty was encountered during the thermal analysis of the test stand due to the 
small number of input parameters and the large number of "unknowns." In order to better 
understand the scope of the problem, a thermal analysis flowchart was developed. Figure 
5-3 shows the test stand thermal analysis flowchart. Data collected from the test stand are 
shown as circle and are used as inputs into the thermal model. From this flowchart it can 
be seen that the problem divides itself into the following three sections: (1) thermal 

network calculations; (2) water channel heat removal calculations; and (3) coupling the 
thermal network calculations with the TOPAZ-II Thermionic Transient Analysis Model 
(TIT AM). TITAM is a computer model developed by the University of New Mexico to 
simulate the operation of the entire TOPAZ-II reactor system under steady-state and 
transient conditions (Paramonov and El-Genk, 1994). The basis for the entire analysis is 
the thermal network. 

1. Thermal Network 

The one-dimensional losses through the test stand are modeled using a network 
method. Figure 5-4 shows a simplified representation of the resistance network used for 
thermal modeling the test stand. The network starts at the outside wall of the cylindrical 
emitter and runs outward to the centerline of the water channel. The three major modes 
of heat transfer encountered throughout the network are radiation, conduction, and forced 
convection. Radiation heat transfer was considered a heat transfer mode through all of the 
test stand gaps when the thermal network was initially developed. It is shown later in this 
section that radiation heat transfer plays an insignificant role in the heat transferred 
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Figure 5-3. Test Stand Thermal Analysis Flowchart. 



Figure 5-4. Thermal Model Resistance Network. 
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through the unregulated and regulated helium gaps, however, radiation heat transfer does 
play an important role in the cesium interelectrode gap. Conduction heat transfer occurs 
through all of the metal components and through the gaps in the test stand; including the 
cesium interelectrode gap, the two unregulated helium gaps, and the regulated helium gap. 
Forced convection occurs into the cooling water flowing through the water channel. 

The basic heat transfer equation used for conduction heat transfer through a 
cylindrical wall is found in Wong (1977). Equation 5.1 gives the equation used for the 
metal components of the test stand: 



Q = (Thot - T C oid) (5.1) 

K c 

where: R c = thermal resistance (KAV) 

T hot = hot wall temperature (K) 

Tcoid = cold wall temperature (K) 



The thermal resistance, R c , is calculated for hollow cylinders using Equation 5.2: 



Rc = 



27tzk 



(5.2) 



where: r 0 = outer wall radius (m) 
r, = inner wall radius (m) 
z = cylinder length (m) 
k = material thermal conductivity (W/mK) 



For the test stand, T hot will be the temperature at the inner cylindrical surface and T cold will 
be the temperature at the outer cylindrical surface of the component of interest. For the 
thermal model calculations, the temperature subscript will indicate the node at which the 
temperature is measured. 
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A cylindrical length of 30.0 cm is used in all of the calculations. This length 
corresponds to the total length of the TISA heater; however, as can be seen in Figure 5-1, 
the length of the TISA heater, the length of the water jacket, and the length of the TFE 
working section are all different. Since there is a 20% increase in cylinder length from the 

TISA heater length to the TFE working section length, the thermal resistance also 
increases by 20%. Although this only results in a few degree Kelvin difference in the AT 
across a single metal component, the cumulative effect of this difference can lead to large 

calculational uncertainty. The shortest of the possible cylinder lengths is used in the 
analysis since it gives a more conservative result. 

The thermal conductivities of all the different test stand components are graphed as 
a function of temperature in Appendix A. Since the thermal conductivity is a function of 
temperature and the metal components have a decreasing temperature distribution from 
the inner wall to the outer wall, an estimate of the thermal conductivity for the entire 
thickness must be made. When entering the thermal model, only one temperature is 
known; therefore, it is not possible to determine the average temperature of the 
component. For the first iteration through the network, the thermal conductivity of the 

component is calculated based on the one known temperature. Because metal is a very 
good conductor of heat, the thermal conductivity is fairly constant over the small AT 
effective across the metal. Based on calculations, using the one known temperature to 

estimate the thermal conductivity results in only a 0.4% difference than if the thermal 
conductivity were recalculated with the lower temperature found after calculating the AT. 
Therefore, less than a three degree Kelvin error is induced into the network and this makes 

recalculation of the thermal conductivity unnecessary. 

The basic equation used for heat transfer for the gas filled gaps is found in El-Oenk 
et al. (1993). This equation accounts for both radiation and conduction heat transfer 
through the gap and is given in Equation 5.3; 
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(5.3) 




where R. c = radiation/conduction thermal resistance (KAV) 



The radiation/conduction thermal resistance, R^, is calculated using Equation 5.4: 



The thermal conductivity for helium as a function of temperature is graphed in Appendix 
A. According to Reid et al. (1977), the thermal conductivity of all gases can also be 
affected by the pressure in the gap, although for pressures greater than approximately 1 
torr, the thermal conductivity only increases by about 1% for a pressure change of 760 
torr. Therefore, for the unregulated helium gap pressures of 200-300 torr, the thermal 
conductivity is considered a function of temperature only. This estimation is valid because 
the mean free path for collisions between atoms is much less than the mean free path for 
collision between an atom and the walls; therefore, an atom does not have to travel far 
before transferring its energy to another atom and eventually reaching thermal energies 
The pressure maintained in the regulated helium gap is around the lower end of where this 
approximation is considered valid. According to the experimental data collected, the 
thermal conductivity of helium in the 1-10 torr range is greatly affected by pressure. This 
unusual phenomenon is discussed in detail later in this chapter. 

The effective emissivity can be calculated by assuming that the view factors from 
both surfaces are equal to 1, and the two surfaces areas are also equal by Equation 5.5 
(Welty et al., 1984): 




ln <T7> ,k 

\ J 



(5.4) 



where: k g = gas thermal conductivity (W/mK) 
d = gap width (m) = (r 0 - r,) 

a = Stefan-Boltzmann constant = 5.67x1 0‘ 8 W/m : K 4 
e e ff = effective emissivity 
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Seff “ 



( 5 . 5 ) 




- 1 



where: s, = inner surface emissivity 
e 0 = outer surface emissivity 



Due to the extremely small dimensions and gap size, both of these assumptions are valid. 



The basic equation used for forced convection heat transfer into the cooling water 
is found in Wong (1977). The equation used for calculating the heat transfer into the test 
stand cooling water is based on Newton's Law of Cooling and is given in Equation 5.6: 



Q = 



D (T w»l| T bulk) 

Kf c 



(5.6) 



where: = forced convection thermal resistance (K/W) 

Tw»ii = coo^ n g water channel wall temperature (K) 
Tbuik = bolk cooling water temperature (K) 



The forced convection thermal resistance, R^., is calculated using Equation 5.7: 



ln(fo 

RfC 27tz(r 0 -n)h 



(5.7) 



where: h = convective heat transfer coefficient (W/m : K) 



The convective heat transfer coefficient is a complicated function of the cooling channel 
geometry, the thermo physical properties of the cooling water, and principally of the flow 
characteristics through the channel. This coefficient and the complexities encountered in 
its definition are discussed in detail later in this chapter. 
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An example of the data set used in the thermal network calculations is given in 
Table 5-1. These data are used in the following discussion to explain the thermal model 
calculational procedure. The data collected between 14-18 August, 1994, in support of 
this thesis are tabulated in Appendix B The data presented in Table 5-1 were taken on 
August 17, 1994. As can be seen there, the TFE efficiency is less than 4%. This is mainly 
due to the interelectrode cesium pressure not being at its optimum value for an input 
power of 3000 W. At this power level, an interelectrode cesium pressure of 
approximately 1 torr would increase the TFE's efficiency closer to 5%. However, for 
thermal modeling purposes, optimum TFE performance is not important. Cesium 
optimization will be covered in detail in the companion thesis written by LT J R. Venable 
that will be available after March of 1995. 



Pttsa (W) 


3,412 


Minimum Collector 
Sleeve Temp. (C) 


354.3 


Paz(W) 


3,003 


Maximum Collector 
Sleeve Temp. (C) 


515.6 


Pou, (W) 


109 


Average Collector 
Sleeve Temp. (C) T 6 


453.7 


Q,o* (W) 


2,894 


Cooling Water Flow 
rate (gpm) 


0.69 


Efficiency (%) 


3.63 


Cooling Water Mass 
Flow rate (kg/s) 


0.04 


Cs Pressure 
(torr) 


0.5 


Cooling Water Inlet 
Temp. (C) 


16.6 


Regulated He 
Gap Press (torr) 


10 


Cooling Water Outlet 
Temp. (C) 


32.2/32.1 



Table 5-1. Example Data Set for the Thermal Model. 
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Since the only temperature measured internally to the test stand is the collector 
sleeve temperature, the thermal network is entered at node T 6 . From this node, 
calculations can be made to determine the temperature distribution either inward toward 
the center of the TFE or outward toward the water channel. In order to make these 
temperature calculations using basic heat transfer equations, the amount of heat 
transferred out through the test stand must be calculated. This transferred heat is equal to 
the portion of the TISA input power that is not converted to electricity in the TFE and is 
designated as Q l0SJ . Equation 5.8 shows the definition for Q loss : 



Qloss P az Pout (5-8) 

In a one-dimensional model for the test stand, Q )0M will travel only radially outward 
through the test stand components to be carried away by the cooling water flowing 
through the water jacket. 

Using the node 6 temperature (T 6 ) and the amount of heat lost by the system (Q 1(m ) 
as a starting point, the thermal model is presented in two sections: (1) calculations from 
the collector sleeve to the water channel, and (2) calculations from the collector sleeve to 
the emitter. Detailed discussions about the critical resistances in the network are also be 
presented. The calculations made with the data from Table 5-1 are given in Appendix C. 

a. From the Collector Sleeve to the Water Channel 

The first resistance encountered in the thermal network is a 1.2 mm thick 
stainless steel collector sleeve. The thermal conductivity of the stainless steel is 
determined using Figure A-3. This thermal conductivity is then used to calculate the 
thermal resistance of the stainless steel using Equation 5.2. Using the calculated Q loiJ as 
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the variable Q, Equation 5. 1 can now be solved for T cold , which corresponds to T„ since all 
of the other equation variables have been defined. 

A 0.05 mm unregulated helium gap is the next resistance in the thermal 
network. The thermal conductivity of this continuum helium gas (gap pressure of 
approximately 200 torr) is calculated using Figure A-4. The effective emissivity is then 
calculated using Equation 5.5. This equation requires the stainless steel emissivity which 
is found from Figure A-6 and the copper emissivity which is found from Figure A-7. Both 
of these emissivities are fairly constant over the temperature range of interest. Using 
Equation 5.4 in combination with Equation 5.3, the value for T g can be solved for 
numerically. If the contribution of the radiation heat transfer is ignored through the helium 
gap, and T 8 is solved for based on conduction heat transfer through the gap only, the 
answers for T 8 are the same within four significant digits. For this reason, radiation heat 
transfer through the unregulated and regulated helium gaps can be disregarded and the 
resistance network can be solved based on conduction heat transfer only. 

The next resistance in the network is the copper conductor. Again 
Equations 5.2 and 5.1 are solved to determine the temperature on the outer cylindrical 
surface of the copper conductor, which corresponds to T,. The regulated helium gap is 
the next resistance in the thermal network. 

b. Regulated Helium Gap 

As discussed earlier in this chapter, the pressure of the helium in the 
regulated gap is maintained between 1-10 torr, depending on the input power level, to 
hold the collector sleeve temperature fairly constant. Since the collector sleeve 
temperatures are measured at approximately the same radial location that the liquid metal 
cooling flows through the actual TOPAZ-II core, holding this temperature constant better 
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simulates the actual conditions in the core. The temperature at this radial position in the 
TOPAZ-II core is relatively constant due to the high heat capacity of the liquid metal 
cooling and the high heat rejection capability of the radiator. 

Due to the low gap pressure, the helium is rarefied which indicates that the 
mean free path for collisions between atoms is much larger than the mean free path for 
atom to wall collisions. Paul Agnew of AEA Technology, United Kingdom, currently 
assigned to the TOPAZ-II Project and working in the area of materials science and 
thermionics, has completed a preliminary analysis on the pressure dependence of helium 
thermal conductivity at these low pressures. The initial results of his work indicated that 
helium thermal conductivity should approach its continuum value at a gap pressure of 
approximately 5 torr vice the 1 torr discussed in Reid et al. (1977). However, 
experimental results indicate that the helium thermal conductivity should be around 30% 
of its continuum value for a gap pressure of 10 torr. Detailed Monte Carlo computer 
programming and investigation of gas atom scattering mechanisms will be required to 
determine a more accurate model of the regulated gap. (Agnew, 1994) 

Although a more exact solution is not possible without extensive work, the 
regulated gap helium thermal conductivity can be better approximated using the empirical 
data collected on the test stand. Since the regulated helium gap is separated from the 
cooling water channel by only 2.0 mm of stainless steel, the temperature of the inner wall 
of the cooling channel is directly dependent on the regulated helium gap thermal 
conductivity. Using the helium continuum thermal conductivity value calculated to be 
0.276 W/mK for this example, gives a cooling channel inner wall temperature of around 
325 °C. The cooling water is pressurized on average to a pressure not greater than 70 
psig Since the cooling system provides water to many different components, the actual 
pressure is dependent on the cooling system load conditions. This cooling water pressure 
gives a initiation boiling point temperature of approximately 140 °C. The wall temperature 
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calculated using the helium continuum thermal conductivity is around 185 °C greater than 
this boiling point temperature. Based on an cooling water outlet temperature of 32 °C, it 
is assumed that the water does not boil in the cooling channel. Figure 5-5 shows a graph 

of the cooling water channel inner wall temperature dependence on the regulated helium 
gap thermal conductivity. If the wall temperature is bracketed around 80 ± 30 °C, the 
regulated helium gap thermal conductivity ranges between 0.073-0.087 W/mK. This 

range is around 30% of the calculated continuum value. 



This calculated range of helium gap thermal conductivities can be used 
along with Equation 5.1 to calculated a range of node 10 temperatures. Using these node 
10 temperatures, the 2.0 mm thick stainless steel water jacket thermal conductivity can be 
determined using Figure A-3. The temperature at the inner wall of the cooling water 




Figure 5-5. Cooling Water Channel Inner Wall Temperature Dependence on 
the Regulated Helium Gap Thermal Conductivity. 
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channel, corresponding to the node 1 1 temperature, is then calculated using Equations 5.2 
and 5.1. These calculated wall temperatures range from 64-122 °C. 

Using an estimated convective heat transfer coefficient and an estimated 
heat transfer area (the method used to estimate these parameters is discussed later in this 
chapter), the bulk water temperature can be estimated using Equations 5.7 and 5.6. These 
calculations yield a range of bulk water temperatures of 6-92 °C. Using an average of the 
measured inlet and outlet cooling water temperatures gives a bulk temperature of 
approximately 24 °C. 

c. Front the Collector Sleeve to the Emitter 

The calculations for the inward portion of the thermal network are similar 
to the calculations made for the outward portion. The first resistance encountered is the 
inner portion of the stainless steel collector sleeve. The same thermal conductivity value is 
used for the collector sleeve as was used in the outward portion. Using Equations 5.2 and 
5. 1, the node 5 temperature can be found. This corresponds to T 5 . 

The next resistance in the network is another 0.05 mm unregulated helium , 
gap. Radiation heat transfer is ignored and the heat transfer calculations using Equations 
5.2 and 5.1 provide the temperature at node 4. Node 4 corresponds to the outer wall of 
the iayer of alumina insulation on the collector. The thermal conductivity of alumina is 
found using Figure A- land the same conduction heat transfer calculations are completed 
to calculate the node 3 temperature, T 3 . This temperature is the interface temperature 
between the collector and the alumina insulation. 

The temperature at node 2, T 2 , can be similarly calculated using the 
thermal conductivity of polycrystalline molybdenum found from Figure A-5. This 
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temperature corresponds to the inner surface of the collector and is an extremely 
important measurement to determine thermionic performance. The cesium interelectrode 
gap is the next resistance in the network. 

iL Cesium Interelectrode (lap 

The cesium interelctrode gap heat transfer mechanism is another critical 
resistance in the network that cannot be easily calculated with only basic heat transfer 
equations. Conduction and radiation heat transfer both play an important part in the 
overall heat transfer through the interelectrode gap. The thermal conductivity of the 
cesium in the gap depends on the gases vapor pressure, the emitter surface temperature, 
the collector surface temperature, and the gap size. (Paramonov and El-Genk, 1994) 

Radiation heat transfer cannot be ignored through the interelectrode gap 
since the temperature of the emitter is high enough to make this mode of heat transfer 
significant (on the order of 1800 K). The effective emissivity for the gap cannot be 
calculated as presented for the helium gaps, but must be determined empirically. A 
detailed discussion on thermal modeling of the cesium interelectrode gap and coupling the 
thermal model with a thermionic emission model can be found in Paramonov and El-Genk 
(1994). 



2. Water Channel Calculations 

In order to calculate the thermal resistance using Equation 5.7, the convective heat 
transfer coefficient is required. Calculating the heat transfer coefficient for the Water 
channel is complicated by the flow channel geometry and by the low water flow rates. 
Figure 5-6 shows the cooling water jacket details including the flow guidance coil. This 
flow guidance coil is a stainless steel tube with an outer diameter of 10.0 mm that fits 
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Flow Guidance Tube in the Water Jacket 



Cooli ng 
Water 




Inlet 



[Flow Guidance Tube 




Flow Channel 




Figure 5-6. Cooling Water Jacket Details. 
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inside of the water jacket to direct the flow in a spiral path to facilitate forced convection 
heat transfer. The cooling water flows around the outside of the flow guidance coil with 
no flow through the coil itself. 

As shown in Figure 5-6, the flow channel cross-sectional area is calculated by 
determining the cross-sectional area of a rectangle with dimensions of 35x10 mm, and 
subtracting off the cross-sectional area of the flow guidance coil (diameter of 10 mm) 
This gives a cross-sectional channel flow area of 2.7x10"* m : . The equivalent hydraulic 
diameter is calculated using the definition from Wong (1977), to be 0.01 1 m. 

Following the definitions and equations presented by El-Genk and Rao (1986), the 
Reynolds number is calculated using this geometry and the given example flow rate to be 
1480. This Reynolds number indicates that the flow falls into the combined laminar flow 
regime. The convective heat transfer coefficient can be defined from El-Genk and Rao 
(1986) and is given in Equation 5.9: 



_ Nu 

D e ( 5 - 9 ) 



where: Nu = Nusselt number 

k = water thermal conductivity (W/mK) 
D e = equivalent diameter (m) 



In the combined laminar regime, the Nusselt number is dependent on the thermophysical 
properties (defined by the Graetz number) and the buoyant forces (defined by the Rayleigh 
number). This flow regime is immensely difficult to analyze with existing correlations. 



Although the convective heat transfer coefficient cannot be determined exactly, 
several existing forced convection correlations can provide a range of possible 
coefficients. The upper limit is determined using the Dittus-Boelter correlation for fully 
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developed turbulent flow in a channel. This gives a convective heat transfer coefficient of 
1012 W/m : K. The lower limit is found by the Seider-Tate correlation for developing 
laminar flow. This correlation gives a convective heat transfer coefficient of 528 W/m : K. 
Therefore, the actual coefficient will fall between the range of 528-1012 W/m : K and these 
limits can then be used in the resistance network to estimate the drop between the cooling 
channel wall temperature and the bulk water temperature. 

Using the cooling water inlet and outlet temperatures along with the measured 
cooling water flow rate, the amount of heat carried away by the water can be calculated 
using Equation 5.10: 



Q m Cp (Tout Tj n ) (5.10) 

where: m = cooling water mass flow rate (kg/s) 
c p = specific heat capacity of water 
T out = cooling water outlet temperature (K) 

T m = cooling water inlet temperature (K) 

The heat carried away by the water, Q, should be equal to the heat lost by the TFE, Q loss , if 
all of the excess heat traveled one-dimensionally out of the test section to be carried away 
by the cooling water. For the example data, the heat carried away by the cooling water is 
approximately 2840 W. This is less than 2% different from the 2894 W lost by the TFE. 
As can be seen from the experimental data presented in Appendix B, most of the data 
collected for the heat removed by the cooling water falls within 10% of the heat lost by 
the TFE. The percent error tends to increase as the input power level decreases. This is 
an expected result since the difference between the inlet and outlet cooling water 
temperature decreases as the power level decreases causing an increase in the expected 
error. 
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By setting Equation 5.10 equal to Equation 5.6, the convective heat transfer 
coefficient can be calculated based on a heat transfer surface area Once again difficulty 
arises in determining the exact heat transfer surface area From Figure 5-6 it can be seen 
that even though no cooling water flows through the flow guidance coil, the stainless steel 
coil is in thermal contact with the stainless steel water jacket The maximum possible heat 
transfer surface area can be determined by summing the inner water jacket cylindrical 
surface area, the outer water jacket cylindrical surface area, and the total surface area of 
the flow guidance coil. The minimum possible heat transfer surface area is equal to the 
cylindrical surface area of the inner water jacket. These calculations are shown in 
Appendix D and give a range of heat transfer area between 0.067-0.203 m : . Figure 5-7 
shows a graph of the heat transfer coefficient plotted against the heat transfer area. Since 




Figure 5-7. Cooling Water Heat Transfer Coefficient Dependence on the 
Heat Transfer Area. 
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the cooling water bulk temperature is known, the cooling channel wall temperature is 
varied from 45 °C to 95 °C. It can be seen from the figure that a wall temperature of 
greater than 95 °C would fall outside of the heat transfer coefficient/heat transfer area 
bands. This information can be used along with Figure 5-5 to further limit the regulated 
helium gap thermal conductivity used in the thermal resistance network. 

3. Thermal Network Coupling With TITAM 

The TITAM computer code provides a means of simulating the operation of a 
TOPAZ-II TFE in the reactor core under both steady-state and transient conditions 
(Paramonov and El-Genk, 1994). Although the test stand is constructed differently than 
the actual core and the thermal and electrical coupling of the TFE with the actual reactor 
systems does not exist, the location of the collector sleeve thermocouples in the test stand 
does provide a means of entering TITAM and having TITAM determine the operating 
characteristics of the TFE. 

The required inputs into the electrically heated TFE component operation portion 
of TITAM are thermal input power, inlet NaK coolant temperature, and NaK coolant flow 
rate. The thermal input power to the TFE is equal to the active zone power, P K , used in 
the thermal network. The average collector sleeve temperature, T 6 , may be used as the 
NaK coolant inlet temperature, and a typical NaK flow rate of 1.3 kg/s is used as the last 
required input to TITAM. From the inputted data, TITAM outputs the radial temperature 
distribution through a TFE in the reactor core and provides the output operating 
characteristics. The TITAM program was executed using the inputs from Table 5- 1 . The 
steady-state output and the relevant data from the thermal analysis are presented in Table 
5-2. 
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Parameter 


TITAM Output 


Thermal Analysis 
Results 


Thermal Power to 
TFE (W) 


3,003 


3,003 


TFE Load Voltage 
(V) 


0.81 


1.38 


Load Current 
(A) 


184.38 


78.89 


Load Electric Power 
(W) 


149.58 


108.87 


TFE Efficeincy 
(%) 


4.98 


3.63 


Cesium Reservior 
Temperature (K) 


576 


528 


Emitter Surface 
Temperature (K) 


1,835.57 


- 


Collector Surface 
Temperature (K) 


765.71 


763.93 


Insulator 
Temperature (K) 


763.75 


761.88 


Coolant Average 
Temperature (K) 


727.58 


726.83 



Table 5-2. Comparison of TIT AM output to Thermal Analysis Results 



The discrepancy in the TFE output power could be the result of the differences 
between operating a TFE in the test stand and simulating the operation of a TFE in the 
reactor core. Another possible source of error are the initial and operating conditions that 
TIT AM uses in its calculations being different from the actual test stand conditions. 
TITAM has a user friendly interface through two input files which can be changed to 
reflect various initial and operating conditions in the reactor core. This includes the ability 
to vary the steady-state interelectrode gap cesium pressure. 
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TITAM initially uses a steady-state cesium pressure of approximately 2.0 torr 
(based on a cesium reservoir temperature of 576 K). The input files were not altered prior 
to running the example data set through TITAM to obtain the above results. When 
TIT AMs initial steady- state cesium pressure was changed to reflect the actual operating 
pressure in the test stand of 0.5 torr (based on a cesium reservoir temperature of 528 K), 
the TFE output power increased from 151 W to 188 W. Although TITAM is being used 
for a purpose that it was not designed for, this trend is contrary to the experimental data 
collected. 

The TITAM temperature distribution correlates well with the thermal analysis 
mainly due to these calculations being based on simple heat transfer principles. This 
indicates that TITAM could be used to estimate the emitter surface temperature based on 
test stand thermal analysis data. 

In order to understand the discrepancies encountered by trying to couple the 
experimental test stand data with TITAM, a detailed examination of how TITAM actually 
performs it calculations and the assumptions under which it operates is required. 
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VI. CONCLUSIONS 



The TOPAZ-II TFE test stand provides an effective means to investigate the 
performance of an individual TFE without the interfering effects and complexity of the 
integrated TOPAZ-II system. A detailed thermal analysis is completed to better predict 
the behavior of a single-cell TFE and to understand the thermal losses associated with the 
test stand components. 

Prior to developing the thermal model, extensive research of the test stand 
construction and materials was necessary. Due to the large volume of Russian 
documentation and the difficulty in translating highly technical data and engineering 
drawings, the development of test stand cross-sectional drawings and the determination of 
component materials and material properties were crucial to the thermal analysis. 

A thermal resistance network has been utilized to determine the thermal loss 
mechanisms through each of the test stand components. Conductive heat transfer through 
the metal components accounts for a very small portion of the overall temperature drop 
from the emitter to the water channel. For the data set used in the example calculations, 
only 7% of the estimated temperature drop from the inner surface of the collector sleeve 
to the bulk water temperature occurred in the metal components. 

Both radiation and conduction heat transfer occur through the regulated and 
unregulated helium gaps, but calculations indicate that radiation heat transfer does not 
contribute significantly to the total heat transferred through these gaps. The pressure in 

the regulated helium gap is varied during power changes to maintain the collector sleeve 
temperature relatively constant. Due to the very low pressure maintained in this gap (< 10 
torr), the helium is rarefied and the thermal conductivity of the gas is not only dependent 
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on gas temperature but is also a function of gas pressure. Estimations based on the 
assumption that the cooling water is not boiling in the channel give a ranse of regulated 
helium thermal conductivities between 0.073-0.087 W/mK at a gap pressure of 10 torr. 
This range is approximately 30% of the thermal conductivity value calculated assuming a 
temperature dependence only. The regulated helium gap is one of the most critical 
resistances in the thermal network. 

The primary mode of heat transfer into the water channel is forced convection heat 
transfer. Determining the amount of heat transferred into the water by forced convection 
requires calculating the convective heat transfer coefficient. This coefficient is dependent 
on the flow characteristics in the channel, the cooling channel geometry and the 
thermophysical properties of the cooling water. The Reynolds number is calculated to be 
1480 which indicates that the flow rate is in the combined laminar flow regime. This is a 
very difficult regime in which to accurately calculate the heat transfer coefficient due to 
the importance of the buoyant, conductive and convective forces in the channel. Based on 
existing correlations, the convective heat transfer coefficient is determined to fall between 
the range of 528-1012 W/m 2 K. 

Another critical resistance in the thermal network is the cesium interelectrode gap. 
Conduction heat transfer occurs through the cesium vapor and can be defined in terms of 
the cesium vapor pressure, the gap size and the surface temperatures of the emitter and 
collector. Due to the high emitter temperatures, radiation heat transfer is significant 
through this gap. The gap heat transfer characteristics are directly affected by the 
thermionic process. Coupling the thermal model with a thermionic model is crucial to 
estimating emitter temperatures based on data collected from the test stand. 

The TITAM computer code provides a means of simulating steady-state and 
transient operation of a TOPAZ-II TFE in the reactor core. TITAM can be utilized to 
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estimate the operating characteristics of a TFE in the test stand by inputting the active 
zone power, average collector sleeve temperature (corresponding to the NaK 
temperature), and an assumed NaK flow rate (based on typical values). The collector 
sleeve temperature calculated by TIT AM is 766 K compared to a temperature of 764 K 
calculated using the thermal network. These temperatures calculations are based on 
simple heat transfer principles which provides good correlation between the two 
calculations. TITAM can also be used to estimate the emitter surface temperature based 
on the input data; however, the thermal resistance of the cesium interelectrode gap must 
be defined before the thermal analysis presented in this thesis can be used to provide an 
estimate of the emitter surface temperature to compare it to. 

Differences in the TITAM output and the thermal analysis results can be attributed 
to several reasons. The most important of these is that TITAM simulates a TFE operating 
in the TOPAZ-II reactor core and not in the test stand. Several other differences include 
cesium reservoir temperature, TFE load voltage, TFE load current, and TFE efficiency. 
To understand the discrepancies between the two results requires a detailed examination 
of the TITAM code and integration of the thermal model developed in this thesis with the 
thermionic model of the code. 
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VII. RECOMMENDATIONS 



The TOPAZ-II TFE test stand provides an excellent vehicle for investigating the 
thermionic characteristics of a single TFE; however, understanding the thermal 
characteristics of the test stand and how its thermal behavior affects the thermionic 
process is crucial in being able to correlate test stand results to the entire reactor core. 
Several improvements to the test stand and additional testing could greatly aid in 
improving the thermal model developed in this thesis. 

A detail investigation of the flow characteristics of the cooling water would 
provide a better estimation of the convective heat transfer coefficient. This could include 
moving the inlet and outlet cooling water thermocouples inside the vacuum chamber 
closer to the test sections water jacket. Using more accurate temperature measurements 
would allow increasing the water flow rate into a better understood and documented flow 
regime. Testing of the cooling water characteristics could be accomplished using a known 
power input into the degassing heaters and would not require the operation of the TFE. 
Thermocouple instrumentation inside the water channel itself, would also prove useful. 

Further understanding the pressure dependence of the regulated helium gap's 
thermal conductivity is critical to the accuracy of the thermal model. Once the thermal 
characteristics of the regulated helium gap and the cooling water channel are well defined, 
a more accurate heat balance can be performed between the collector sleeve and the 
cooling water. 

Exploring the end effects at the upper and lower helium chambers would provide a 
better understanding of the 2-dimensionality of the thermal system. Experiments could 
easily be performed without the operation of the TFE to examine how these end chambers 
are affecting the axial temperature distribution along the collector sleeve. Operating the 
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heaters for one or both of these chambers during operation of the TFE provides an 
additional heat input into the thermal model; therefore, the input power to these heaters 
should be recorded. 

Experimental information is only valuable if the accuracies of the instrumentation 
and data acquisition system are completely defined. The integration of Russian and 
American equipment is inherent to the TOPAZ project and increases the difficulty in 
developing an accurate model for the test stand. Improved documentation of the test 
stand measurement errors should accompany any future experimental work. 

Experimental data should be collected from the test stand at a predetermined TFE 
output voltage. This voltage level should be based on supplying a specified bus voltage 
(i.e. 28 VDC or 30 VDC). This must be a consideration since operating the TFE at the 
maximum output power (or maximum efficiency) might not produce useful voltage to 
supply a power distribution bus. 

TITAM is a valuable tool that should be modified further for use with the test 
stand. Utilizing the thermionic emission model from TITAM along with the thermal 
model developed in this thesis can provide a detailed understanding of the thermal 
characteristics of the test stand. This would allow for a more direct correlation of the test 
stand results with the operation of the overall TOPAZ-II reactor system. 
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APPENDIX A. THERMAL PROPERTIES 



The thermal conductivities and emissivities of various test stand components are 
graphed and are defined in terms of temperature. Most of the information was available in 
Paramonov and El-Genk, 1994, but some was derived from reference book data tables. 
References are provided in the titles of each figure. A polynomial curve fit of the 
tabulated data accomplished using Mat head Plus 5.0 (MathSoft, Inc., 1994). 
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Temperature (K) 

Figure A-l . Thermal Conductivity of Alumina (A1 2 0 3 ). 
(Paranomov and El-Genk, 1 994) 
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Figure A-2. Thermal Conductivity of 99.99% pure Copper. 
(Powell et al., 1966) 



Figure A-3. Thermal Conductivity of AISI-321 Stainless Steel. 
(Peckner and Bernstein, 1 977) 
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Figure A-4. Thermal Conductivity of a Helium Continum. 
(Andrews and Biblarz, 1981) 



(Paraiiomov and El-Genk, 1 994) 
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Figure A-6. Emissity of AIS 1-3 21 Stainless Steel. 

(Peckner and Bernstein, 1977) 
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Figure A-7. Emissity of 99.99% pure Copper. 

(Baumeister, 1978) 
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Figure A-8. Emissity of Alumina (A1 2 0 3 ). 
(Baumeister, 1978) 



APPENDIX B. EXPERIMENTAL DATA TABLES 



Data were collected in support of this thesis during two separate experimental 
sessions The first occurred between 15-26 March, 1994, when data were collected over a 
TISA heater input power range of 1000 W to 3600 W. The primary purpose of the runs 
completed during this period was to establish the operating characteristics of the test stand 
and its associated TFE from normal operating range to very low input power levels. This 
information is being used in a follow-on thesis to examine the thermionic characteristic of 
a TFE in the test stand. 

The second experimental session occurred between 14-18 August, 1994 Prior to 
this session, the entire test stand data acquisition system was converted from a computer 
hardware and software system that was furnished by the Russians to a more contemporary 
system that was designed and assembled by the US, TIP personnel. This system update 
provides for much easier data acquisition and test stand control. Although the data 
acquisition system changed between the two experimental procedures, the test stand 
operating conditions were maintained as close as possible to allow for data correlation 
between the two collection periods. The primary purpose of the second experimental 
procedure was to support thermal model calculations by investigating the relationship 
between regulated helium gap pressure and collector sleeve temperature and to collect 
data a various power levels for an optimized cooling water flow rate. Relevant data to the 
thermal analysis are tabulated in this appendix. 

The appendix is separated into three time periods, each time period consisting of 
three tables. The shaded columns contain data that were recorded from the test stand and 
the unshaded columns contain derived data. All of the data are collected at a cesium 
reservoir temperature of approximately 255 °C (around 0.5 torr of pressure). 
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The second column of the first table is TISA active zone power and is determined 
by multiplying the TISA input power by 88%. The fifth column is TFE output power and 
is determined by multiplying the TFE output voltage by its output current. The column 
indicated as column "A" is the power loss from the TFE. These values are found by 
subtracting the TFE output power from the TISA active zone power. TFE efficiency in 
column seven is the result of dividing the TFE output power by the TISA active zone 
power and multiplying by 100. 

Column eight of the first table is the average collector sleeve temperature and is 
the same value as in column five of the third table. Performing a heat balance on the 
cooling water provides the data in the columns indicated as column "B" and "C". There is 
one column for each of the two output thermocouples. These values are the same values 
as given in columns 12 and 13 of the third table. The % error in columns 1 1 and 13 of the 
first table are the percent difference between the cooling water heat balance (columns "B" 
and "C") and the power balance (column "A") 

The second table presents the 11 collector sleeve thermocouple measurements. 
See Figure 5-1 for the relative axial location of the thermocouples along the collector 
sleeve. Columns 2 through 6 of the third table provide relevant information concerning 
the thermocouple measurements. Column five gives the average of all 1 1 collector sleeve 
thermocouples and column six gives the average of thermocouples T57 through T66. This 
excludes the end point temperature measurements of T56 and T67. 

Column eight of the third table contains the measured water flow rate (in gpm) 
converted to a cooling water mass flow rate (in kg/s). Columns 12 and 13 of the third 
table indicate the results of the cooling water heat balance (utilizing Equation 5.10). 
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The tables dated 13-14 August 1994 are taken over a TISA input power of 
1000-3600 W. At each power level the regulated helium gap pressure was adjusted to its 
optimum value. This maintains a relatively constant collector sleeve temperature over the 
entire power range. The high % error for the 1000 W data run is due to the small 
temperature difference between the cooling water inlet and outlet. This indicates how 
sensitive the cooling water heat balance is to this temperature difference. 

The tables dated 14-16 August 1994 provides data collected over three power 
ranges with varying regulated helium gap pressure at each power level. The first table 
shows that for a constant input power level, the collector sleeve temperature decreases as 
the regulated helium gap pressure increases. This proves how important understanding 
the pressure dependence of the regulated helium gap thermal conductivity is to the thermal 
analysis of the test stand. 

The tables dated 17-18 August 1994 provides a set of data runs under similar 
conditions as presented in the first set of tables. A complete set of all data acquired during 
both experimental periods is archived at the TIP facility in Albuquerque, New Mexico. 
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Experimental Test Stand Data taken 14-16 August 1994 
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APPENDIX C. THERMAL NETWORK CALCULATIONS 



The following calculations and the calculations in Appendix D were completed using 
Mathcad Plus 5.0 (Mathsoft, Inc. 1994). This calculational software uses a colon followed 
by an equal sign (:=) to mean "defined by" vice an equal sign (=) to mean "equal to". All 
variables must be defined prior to being used in an equation. After a variable is defined, an 
equal sign can be used to determine the variables value. 

A. THERMAL MODEL INPUTS 

1. INPUT POWER TO THE TISA HEATER 

Q TISA ; 34 12- watt 

2. POWER TO THE ACTIVE ZONE 

Q az = Q TISA 0-88 Q az = 3.003 -10 3 -watt 

3. TFE OUTPUT POWER 
P out - 108.9 watt 

4. POWER LOSS FROM TFE 

Q loss = Q az " ^ out Qloss - 2.894*10 -watt 

5. RADIAL DIMENSIONS FOR THE TEST STAND 



= 0.0098 m 


r 4 


^ 0.0119 m 


r 7 


= 0.01495 m 


r 10 


= 0 0300 m 


= 0.0103 m 


r 5 


= 0.01195 m 


r 8 


= 0.0150 m 


r ll 


= 0.0325 m 


= 0.01175 m 


r 6 


-- 0.01370 m 


r 9 


= 0.0295 m 


r 12 


= 0.0375 m 



6. LENGTH OF THE TEST SECTION 

z = 0.3 m 

7. DIMENSIONAL CONSTANT FOR CALCULATIONS 

c = 2-Tt-z c = 1.885-m 
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8. COLLECTOR SLEEVE AXIAL TEMPERATURE DISTRIBUTION 

T 56 = 424.49 K T 59 = 483. 9 K T 62 433.73 K T 66 = 362.95 K 

T 57 = 515.64 K T 60 = 482. 57 K T 63 = 439.32 K T 67 = 354.29 K 

T 58 = 498.46 K T 61 = 490. 12 K T 65 = 376.42 K 

9. AVERAGE COLLECTOR SLEEVE TEMPERATURE 

The average collector sleeve temperature is calculated without using the end point 
temperatures, T(56) and T(67). 




T 6 •= T 6c * 273.15-K T 6 = 726.829-K 



B. FROM THE COLLECTOR SLEEVE OUTWARD 



1. STAINLESS STEEL COLLECTOR SLEEVE 



a. Thermal Conductivity of Stainless Steel 




b. Temperature Drop Across Collector Sleeve 




AT ss = 6.957-K 



c. Temperature at node 7 



T 7 =T 6 -AT SS T 7 =719.872-K 
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2. UNREGULATED HELIUM GAP 



a. Thermal Conductivity of a Helium Continuum 



k He2 = 5 - 5 ' 10 



4 / watt 


T _ _ 


9.13 10 7 -l 

\ 


watt ^ 2 




\m-K 2 


) 7 


_ k 3 i 7 

m K. / 




/ watt \ 


T 3 . 


1 264 10' 13 


/ watt \ 4 




' m K 4 / 


1 7 


nvK 51 7 




/ watt 
m-icj 


" k Hel 


T k He2 


k He = 0.282- 


watt 

m-K 



b. Temperature Drop Across the Unregulated Helium Gap 
r 8 



AT 



He 



Q loss ^ nl 

VtI 



c k 



He 



AT He = 18.147-K 



c. Temperature at Node 8 



T 8 = T 7 - AT He T g = 701.725-K 



3. COPPER CONDUCTOR 



a. Thermal Conductivity of Copper 



k Cu = 422.286 — a - - 6.71429 10 



m-K 



watt \ 

\m-K 2 j 



T g k Cu =375.17 



watt 

m-K 
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b. Temperature Drop Across the Copper Conductor 



AT 



I r, 

Q loss' 1 " “I 
r 8 



Cu 



ck Cu 



AT Cu =2 767 ’ K 



c. Temperature at Node 9 

T 9 --T 8 -AT Cu T 9 = 698.958*K 

4. REGULATED HELIUM GAP 



The thermal conductivity of the helium in the regulated gap is first estimated using 
the equation for a continuum gas. 



a. Thermal Conductivity of a Helium Continuum 



k He3 1 9.376- 10 | 



watt \ 
m-K 2 j 



•T 9 - 9.13-10 



-7 / watt \ 2 

3 T 9 

y m-K 3 / 



k He4 - 5.5-10 



-10 



watt \ 
\in-K 4 ) 



•T 9 - 1.264-10 



. i 3 / watt \ 4 

T 9 



,m-K' 



. ~ . , -2 /watt\ , , . ^ watt 

k Hereg = " 9 06 10 ' ~ ~ ^ k He3 " k He4 k Hereg “ a276 *“ T7 



b. Temperature Drop Across the Regulated Helium Gap 



AT 



Hereg 



Q loss ln 


io\ 




l r 9 / 



ck Hereg 



AT Hereg = 93366 * K 
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c. Temperature at Node 10 Based on a Helium Continuum in the Gap 



T lOcont = T9 ^ Hereg T lOcont 605.592-K 

This temperature will lead to a water channel inner wall temperature of around 
325 °C. Since this is approximately 200 °C higher than the boiling point of water, the 
helium thermal conductivity must be much lower than its continuum value. From Figure 
5-5, the helium thermal conductivity is estimated to be 0 08 W/mK for a wall temperature 
of 80 °C. Using a temperature range of 80 + 30 °C gives a low thermal conductivity of 
0.073 W/mK and a high thermal conductivity of 0.087 W/mK. The temperature drop 
across the regulated helium gap is recalculated based on the low, estimated, and high 
thermal conductivities. 



d. Estimated Temperature Drop Across the Regulated Helium Gap 



. „ watt 

^ Helow = — — ^T 



Q loss * n 



10 



m-K Helow 



c k 



Helow 



Q loss'* 11 



10 



watt 



k Heest = ® ATpj ee$t 



m-K 



ck 
Q loss * n 



Heest 
r 10' 



watt 



^ Hehigh = 0 087 m .K Hehigh 



ck 



Hehigh 



AT Helow = 353.441 -K 



AT Heest ~ 322.51 5 -K 



AT Hehigh =296.565-K 



e. Estimated Temperature at Node 10 



lOlow = T 9 ’ AT Helow 


T lOlow = 345.517-K 


10 = T9 * AT Heest 


T 10 =376.443 -K 


lOhigh = T 9 AT Heh j gh 


T lOhigh =402.393-K 
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5. Stainless Steel Water Jacket 



a. Thermal Conductivity of Stainless Steel 

^ wjlow = ^^'(nvic) + lOlow ^ wjlow ~ *^^ 



. , „ „ / watt \ n I watt 

k * = 12.0*1 + 0.01 

^ \m-K/ 1--2 



TlO 



wjhigh 



= 12 . 0 - (- a - +■ 0.01 
m-K 



m-K / 
watt 



, ,r~,s, watt 

k »-’ 15 764 m-K 



.m-K 



2 7 lOhigh ^wjhigh 16.024 



b. Temperature Drop Across the Water Jacket 
' r ll 



Q loss'* 0 
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ck 



Q loss'* 0 ! 
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'r 1 1 \ 
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It 
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\ r 10/ 



c k wjhigh 

c. Temperature at Node 11 



AT wjlow =795 * K 



AT vyj = 7.795 -K 



AT wjhigh = 7668 * K 



11 low ' 7 lOlow ^ 7 v^low 


T 1 How = 337.567-k 


II 
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1 
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H 


T jj = 368.649-K 


11 high ' 7 lOhigh " ^ wjhigh 


7 1 lhigh = 394.724-K 



watt 
m K 



watt 

m-K 
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These calculations yield a water channel inner wall temperature of between 
64-122 °C. Using the Dittus-Boelter correlation convective heat transfer coefficient of 
1012 W/m 2 K as the upper limit and the Seider-Tate correlation convective heat transfer 
coefficient of 528 W/m 2 K as the lower limit, the last resistance in the network can be 
calculated and a bulk cooling water temperature can be calculated. 



6. WATER CHANNEL 



a. Convective Heat Transfer Coefficient and Heat Transfer Area 

2 



^ low ~ ^255 
h 



watt 

m 2 K 



htgh » 1012- 2 



watt 



m K 



A htlow s 0 07 m 
A hthigh = 0125 m 

A htlow * A hthigh 
A ht ; - 



A h t = 0.098-nC 



b. Temperature Drop from the Wail to the Bulk Fluid Temperature 

= 56.209-K 



AT 





Q loss 


wlow 


wlow 






A ht' h low 






Q loss 


whigh 


whigh 


A ht'l 1 high 



= 29.327-K 



c. Bulk Water Temperature 



^bulklowl 1 ^ 11 low ' wlow 


^"bulklowl 281.357‘K 


^bulklow2 = ^ 11 low ‘ whigh 


T bulklow2 ~ 308.24-K 


T bulkl = T 11 - AT wlow 


T bulk 1 = 312.439-K 


^ bulk2 = ^ 1 1 ‘ whigh 


T bulk2 = 339.322-K 



^bulkhighl = ^ 1 lhigh “ wlow ^bulkhighl 338.515-K 

^bulkhigh2 = ^ llhigh " A ^whigh ^bulkhigh2 = 365.398-K 

These calculations give a bulk cooling water temperature range of 8-92 °C. 
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C. FROM THE COLLECTOR SLEEVE INWARD 



1. STAINLESS STEEL COLLECTOR SLEEVE 



a. Thermal Conductivity of Stainless Steel 



, , . . / watt \ . . , / watt . _ 

k ss = 12.0-1 I 0.01-1 1-T , 



mK 



iin-K' 



, , ^ watt 

k gg — 19.268 

m-K 



b. Temperature Drop Across Collector Sleeve 



AT 



SSi 



, r 6' 
Q loss* 0 

\ r 5i 



ck SS 

c. Temperature at node 5 



AT ssi = 10.888-K 



T 5 = T 6 + AT 55 } 



Tc = 737.717-K 



2. UNREGULATED HELIUM GAP 



a. Thermal Conductivity of a Helium Continuum 



. , .-4 / watt \ ^ ,«-7 / watt \ ^ 2 

k He li - 9.376-10 • - -T 5 - 9.13-10 • -T 5 



,m-K 



im-K J 



He2i 



- 5.5-10 



- io / watt \ ^ 3 



im-K 



T 5 j - 1.264-10 



13 / watt \ 4 

. )' T 5 
tm-K' 



k He i = -9.06 10 2 (— a £) + kHeli + k He2i k Hei “ 0 288 

\m-K/ m-K 
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b. Temperature Drop Across the Unregulated Helium Gap 



AT 



. , ; r 5 i 

Q loss * n l I 
r 4 



Hei 



c k 



Hei 



AT Hei = 22.382-K 



c. Temperature at Node 8 

T 4 = T 5 + AT Hej T 4 = 760.099-K 



3. ALUMINA INSULATION 



a. Thermal Conductivity of Stainless Steel 



watt 



watt 



k aluml ' 74 9 -r: - 0.178 r T 4 ♦ 1.75 10 



- 4 watt 



m K 



m-K 



m-K 



3 



, nno 1A -8 watt 3 , ,,, , n -n watt 4 

k alum2 = _7 - 7810 7 t 4 +l311 ° 7 T 4 

m-K 4 m-K 5 



k alum r k aluml + k alum2 



, , . _ , . watt 

k alum = 10916 ~ 

m-K 



b. Temperature Drop Across the Alumina Insulation 

I 

Q loss 4n 



r 4 \ 



AT 



alum 



\ r 3/ 



c k alum 



AT alum = 1784 - K 



c. Temperature at node 3 

T 3 ■ T 4 + ^ alum T 3 = 761.883 -K 
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4. COLLECTOR 



a. Thermal Conductivity of Polycrystalline Molybdenum 



k polyMo • ' ' 510 ‘ 2 T7T T 3 k polyMo= 98 572 ' 



watt 
m K 



b. Temperature Drop Across the Collector 



Qloss * n r 

AT coll =-7—^ ATcoH-Z.OSl.K 
c K polyMo 

c. Temperature at Node 2 

T 2= T 3* AT coll T 2 =763.934-K 
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APPENDIX D. WATER CHANNEL CALCULATIONS 



A. WATER CHANNEL HEAT TRANSFER AREA 

1. Dimensional Constants 



r|| = 0.0325 m Inner Cylinder Radius in m 

r J 3 = 0.0425 m Outer Cylinder Radius in m 

z = 0.326 m Water Jacket Length in m 

r c - 0.005 m Radius of the Flow Guidance Coil 

d c j = 0.075 m Test Stand Diameter to center of the Water 

Channel in m 

9 - 13 deg Cant Angle of Flow Guidance Coil in degrees 

N = 6.5 Number of Revolutions of the Flow Guidance 

Coil 



2. Inner Cylinder Surface Area in m 2 : 

A j - 2-7CT J J z 
A j = 0.067* m 2 

3. Outer Cylinder Surface Area in m 2 : 

A o = 2k ' t 13 z 
A 0 = 0.087*m 2 

4. Surface Area of the Flow Guidance Coil in m 2 : 



d e £* = Effective Test Stand Diameter based on the Flow 

cos (9) Guidance Coil Cant Angle in m 

d e ff = 0.077*m 

L - N-7t -d e g* Length of the Flow Guidance Coil in m 

L = 1.572*m 

C c - 2-7tr c Circumference of the Flow Guidance Coil in m 

C c = 0.03 l*m 

A c .= L C c Surface Area of the Flow Guidance Coil in m 2 



A c = 0.049*m 2 
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